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Abstract: We investigate heavy-quark production as a function of the rapidity interval 
between two heavy quarks in hadronic colhsions. We compare the results relevant to 
bottom production at the Tevatron and at LHC, obtained using exact leading-order and 
NLO pQCD production, as well as the contribution of the 46 channel with and without the 
addition of BFKL gluon radiation. 



Keywords: 



Standard Model, QCD, Heavy Quarks, Hadronic Colliders . 



1. Introduction 



One of the most important processes at high-energy hadron-hadron colhders is the pro- 
duction of heavy quarks. Bottom and top quark production, for example, provide not only 
many tests of perturbative QCD, but also some of the most important backgrounds to 
new physics processes. Not surprisingly, therefore, such heavy-quark production has been 
extensively studied in the literature (see e.g. Ref. for a review) and the phenomenology 
at the Tevatron and the LHC has been evaluated in great detail. 

In the kinematic region in which the transverse momentum of the heavy quark Q is 
of the same order as its mass mg, the leading-order contribution to the inclusive heavy- 
quark production cross section comes from the partonic subprocesses in which a QQ pair 
is produced, gg, qq — > QQ. The next-to-leading order (NLO) corrections to these processes 
have been available for quite some time now |2|, |^, ^, ^. They are numerically important, 
particularly for h quarks, where they can result m. a K factor as large as two. 

At the parton level, these large radiative corrections to the total rates are easily iden- 
tified as coming from production near threshold, s ~ 4mQ {s being the partonic centre- 
of-mass energy squared). When folding partonic cross sections with parton distribution 
functions (pdfs) to get the observable rates, the threshold region is especially relevant in 
those cases in which the total hadronic energy ^/S is of the same order as the quark mass, 
as for example for top production at the Tevatron, or h production at fixed-target facilities. 
Potentially large logarithms appear in the perturbative expansion, and these need to be 
resummed to all orders. In practice, however, this resummation only marginally increases 
the NLO predictions (see e.g. Ref. |^]). 

Total partonic rates can also receive large contributions from the high-energy region 
s ^> 4mQ, complementary to the threshold region. As discussed in Ref. |2[, this is due to 
those partonic subprocesses that feature a gluon exchange in the t-channel; this happens 
for gg QQg and qg QQq, and it is peculiar to the NLO computations of quark pair 
production, as opposed to Born-level predictions, in which only fermions are exchanged in 
the t channel. It must be stressed that at the hadron level this enhancement is diluted by 
the fall-off of the pdfs at large x values @, |^ . 

A gluon exchange in the t-channel is also present at 0(a|) in the reaction gg QQQQ, 
which is the Born-level contribution to this four-quark process. This is interesting, since the 
t-channel gluon exchange leads to properties fairly similar to those relevant to the Mueller- 
Navelet dijet cross section which is used to study the high-energy limit of QCD in 
which the energy dependence of the lowest-order cross section is enhanced by BFKL-type 



logarithmic corrections ||10|, |ll|, |12 



The dominance of the gluon exchange in the t-channel implies that the 4(5 channel 
is perturbatively suppressed only by a factor of ag with respect to pair production at 
high energies. Although this still prevents us from a straightforward use of AQ production 
to detect BFKL signals, we can, however, observe that in the high-energy regime the 
kinematics of the 2Q and 4Q production channels are rather different. The former is 
dominated by those configurations in which the QQ pair recoils with large rapidity against a 
fast light parton. On the other hand, the AQ system will predominantly be produced in two 
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QQ pairs, rapidly moving away from each other; the relative rapidity of each pair is small 
compared to the separation in rapidity of the two pairs. Therefore by selecting particular 
kinematic configurations it may be possible to relatively enhance the AQ contribution and 
find signatures of BFKL. This is the main focus of our study. 

In order to define a proper set of observables, we require for each event to tag (at least) 
two heavy flavours (in any possible combination: QQ, QQ, or QQ), which we denote by 
Qi and Q2, separated by a large rapidity interval, Ay = |yQ^ — yQ^\ » 1. In this way, we 
should cut off the configurations that dominate pair production in the high-energy regime, 
without losing too many events in the AQ channel. We aim at studying whether this is 
the case or not, specifically in the regions accessible to the detectors at present and future 
colliders, by comparing the predictions for 2Q and AQ production processes. We stress 
that our set of observables is based on a double Q tagging, which in fact is already used to 
study QQ correlations in heavy-quark pair production. In this paper, we shall not correct 
our results for tagging efficiency. 

The high-energy limit of AQ production can be considered as a reformulation of the 
standard Mueller-Navelet dijet case. What we are doing here, in effect, is replacing each 
Mueller-Navelet forward jet (with pT > -Prmin) by a QQ pair. In fact, by identifying 
the rapidity of each pair with the rapidity of the tagged quark in the pair, we have s = 
Attiq^ cosh^ y*, where w-g^ = + tuq is the squared heavy-quark transverse mass, and 
y* = (Ay)/2. The formula above, relating the large-s to the large-y* region, is customary 
in Mueller-Navelet arguments. The differences between jet and heavy-quark production are 
easy to find: whereas in the dijet case it is -Prmin that regulates the infrared singularities 
at t = 0, here it is the heavy-quark mass mq. The analogue of the P'Xmin^ii ~^ constant 
behaviour of the leading-order dijet cross section at large dijet rapidity separation Ay is 
the rnqCFiQ — > constant behaviour of the AQ heavy-quark cross section. The effect of the 
(leading logarithm) BFKL corrections is the same in both cases: the partonic cross sections 
increase asymptotically as exp(AA7/) where A = 4 log 2Ncas/TT and Ay is either the rapidity 
separation of the dijets in the Mueller-Navelet case, or the rapidity separation of the two 
QQ systems in the present context. 

Another process of potential interest in the high-energy limit is QQ + 1 jet production. 
In this case the partonic subprocesses gg — > QQg and qg — > QQq, which feature a gluon 
exchange in the t-channel, are 0(a|) at the Born level. This can also be considered as a 
reformulation of the standard Mueller-Navelet analysis, where only one of the forward jets 
is replaced by a QQ pair. 

This paper is organized as follows. In Section |2| we compute analytically the high- 
energy limit of the gg QQQQ cross section. In Section |3|, we describe how to include the 
resummation of BFKL logarithms, through Monte Carlo methods. In Section ^, results for 
2Q and AQ channels are compared, at the Tevatron and LHC energies. We also consider 
the case of QQ + 1 jet production. Finally, in Section |5| we present our conclusions. The 
appendices collect some useful formulae. 
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(a) 




Figure 1: Amplitudes for QQ production in gg fusion. Figure (a) represents the leading-order 
term. Figures (b), (c) and (d) are examples of the diagrams that contribute to the NLO term. 
Figures (e) and (f) represent the 4 Q contribution to the NNLO term; figure (e) ((f)) is an example 
of a diagram with quark (gluon) exchange in the t channel. Figure (f) constitutes also the leading 
term of a BFKL gluon ladder, and figure (g) represents the first rung of it. 

2. The high-energy hmit 

In the high-energy limit, the Ay distribution for QQ production can be written schemati- 
cally as 



da. 



QQ 



Ay 



~ a 



j=0 



1 

as 



3=0 



j=Q 



(2.1) 



where L = log(s//i^) ~ Ay is a large logarithm, and the quantity /i^ is a mass scale 
squared, typically of the order of the crossed-channel momentum transfer and/or of the 
heavy-quark masses. The first sum in Eq. (|2.1| ) is a fixed-order expansion in as starting at 
0{al) (the Born processes qq,gg — > QQ), which collects together the contributions that 
do not feature gluon exchange in the crossed channel between the heavy quarks. The aoo 
coefficient is the leading-order term, which for gg fusion is depicted in Fig. |^(a); the agi 
coefficient is the NLO term (specimen diagrams are given in Fig. |^(b-d)). An example of a 
4Q contribution to the ao2 coefficient is given in Fig. ^(e). The aoj coefficients behave like 
1/s, or equivalently exp(— Ay), modulo logarithmic corrections.* In Eq. ( |2.1| ), the second 

*The ao2 coefficient may also contain terms that behave like l/{\/lfiv/) and arise from the interference 
between diagrams with gluon exchange in the crossed channel and diagrams with quark exchange in the 
crossed channel. 
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and third sums collect the contributions which feature only gluon exchange in the crossed 
channel between heavy quarks, the second (third) sum resumming the BFKL (next-to- 
)leading logarithmic corrections. Fig. |l](f) represents the zeroth-order term, and Fig. ||(g) 
contributes to the first-order term, of the second sum. The ay and a2j coefficients behave 
like l/z^w) contrast to the 1/s behaviour of the aoj. The ellipses of Eq. ( |2.1D refer to 
logarithmic corrections beyond the next-to-leading accuracy. Thus, it is clear that the 
second and third sums of Eq. ( |2.1| ) will eventually dominate over the first sum in the 
asymptotic energy region s — > oo. In Sections |2| and ^ we will analyse the second sum of 
Eq. ( |2.l[) in the region s ^ /x^^, by computing the aij coefficients in the high-energy limit. 

Details of the calculation for the production of four heavy quarks, via the sub-processes 
99 ~^ QQQQ and qq —>■ QQQQ, are presented in Appendix In the high-energy limit, 
we require that any two Q's (no distinction between Q and Q is necessary) are produced at 
large rapidity separation. Then the production process is dominated by the sub-processes 
for which the tagged Q's are separated by gluon exchange in the crossed channel. Of the 
above two sub-processes, only gg — > QQQQ features gluon exchange in the crossed channel. 
With the kinematics of the high-energy limit, 

VQi - » VQs - VQi ' PQii^ - PQ2i_ - PQ:ii- - PQii_ ' (2-2) 

the amplitude for gg QQQQ factorises as 

4^2 _ _ 
\^gagt^QiQ2Q3Qf = [l^^^ iPa,PQf ,Pq/, q)!"^"^ {Pb,PQ„PQ,; -q)\ , (2.3) 

where Qf{Qf) and Qb{Qb) are the quarks (anti-quarks) produced forward and backward, 
respectively. In Eq. (^I3|), 

i=q^ = {Pa-PQ,-PQ/ (2.4) 

is the momentum transfer. The impact factor J*^*^ is calculated in Appendix ^ starting 
from the amplitude for — > QQi with i = q,g and using high-energy factorisation. The 
result is given in Eq. ( |B.4| ), summed (averaged) over final (initial) colours and helicities. In 
the kinematics of ( p.2[ ), the exact parton momentum fractions ( [A.(i| ) are well approximated 
by 

^a- ^ , X,- . (2.5) 

Using Eq. ( |2.3D , we can write the cross section for heavy-quark production as 
da 

/ 2 N r 2 N I^^{qa)I^^{qb) (^^(qax - ^b^) ,^ r\ 

= X^ga/AiXa, f^Fa) Xb9b/B{Xb^ f^Fb) ^ 4(A ^4 2 2 ^ ^' i^'^) 

27r4 47r fg^^^ 2 



t Contributions like the one in Fig. |l](c), which feature gluon exchange in the crossed channel but not 
between heavy quarks, are not systematically resummed in Eq. (2.1), and are thus implicitly included in 
the first sum. They contribute, however, to the leading order for QQ + 1 jet production, where a gluon 
is exchanged in the t-channel between the jet and the QQ pair, and they constitute in that case the Born 



term of the BFKL ladder. We will consider QQ + 1 jet production in Section 4.2 
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with momentum transfers Qa = Pa — PQi — Pq^ lb = PQ3 + Pq^ ~ Pb, and where 
9a/A{^a^ ^^'Fa) is the pdf for the gluon ga, and analogously for gi,. We use the notation p_L 
to denote a transverse momentum vector. 

However, in Eq. (2.6) energy and longitudinal momentum are not conserved. The par- 
ton momentum fractions in the high-energy limit, and x^, underestimate the exact ones, 



Xa and Xfe, Eq. (A. 6) and accordingly the values of the pdfs are overestimated. Thus for the 
numerical applications of Section ^ we will use the factorised form ( |2.6| ) of the production 
rate, but with x^ — > Xa and x^ Xf,- This modification is particularly important when 
BFKL evolution is considered. 

The above results must be integrated over the phase space of the final-state particles 
in order to get physical results. In the high-energy limit, the phase space (|A.2D can be 
factorised into the phase spaces for the two impact factors. 



dVA 



n 



X (27r)2 6^{pQ^^ + Pq^^ + pQ3^ + pQ^J , 
where we have used light-cone coordinates = {p^ ibp'^)/\/2. Fixing 



Zb 



PQi+Pq2 

the phase space ( p. 7)) can be rewritten as 

11/ dza d^PQij 



dVA 



Pq,+Pq, 



PqJ 

(2.7) 
(2.8) 



(47r)2 2s V^a(l-^a) (2vr)2 jyzb{l-Zb) (27r)2 
(2vr)^r(qa^ -qb^), 



(27r)2 (27r)2 

with centre-of-mass energy s = 2p^p^ . Note that Eq. (|2.9| ) is written in such a way as to 
be immediately generalizible to the emission of a BFKL gluon ladder between the impact 
factors. 



Using Eqs. (|2.3| ) and (2^) in the expression for the cross section given in ( |A.lJ ), we 
obtain 

d'^Cla^d'^Clb^I{(laJI{(lbJ . 2,2, 



daipaPb^PQ„PQ,PQ,PQ,] - ^^^y ^2vr)2 q2 
where the integrated impact factor is 



(27r)^5^(q,^ -qb^), (2.10) 



dx d^p± 



/««(x,px;qx) 



(2.11) 



4tt Jo x{1 - x) (27r)2' 

with I'^^ given in Eq. ( [B.4| ). The integral is explicitly performed in Section [B.l| , where it 
is expressed in terms of a function g, Eq. (B.14), of the dimensionless ratio ^ = q'j_/mQ. 
Then using Eqs. ( |B.10|) - (|B.14D , the total integrated cross section (|2.10| ) becomes 



93 



— — j — / — — 9\ia)9{ib)Kia - ib) 



(2.12) 
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Figure 2: Partonic cross section for QQ and QQQQ production. The histograms show the exact 
leading-order results, i.e., the exact matrix elements integrated over the exact phase space. The 
diamonds are obtained integrating the high-energy limit of the matrix element, Eq. ( ^.3|) , with the 
exact phase space. The patterned red line is Eq. (2.13), representing the asymptotic limit. For 
comparison, the Born gg, qq — > QQ contributions are also shown. The coupling as is set to one. 
Note that the kinematic limit for fe-quark production at the Tevatron is at F « 10.6. 



Note that even though according to Eq. ( |B.13| ) the function g{^) grows logarithmically 
with ^ as ^ — > oo, the integral in ( ^.12|) is finite and gives |13| 



1 



vrm^ - 1 



23iV2 277 /175C(3) 19 
~81 486 ^ V 576 2887 iVf 



Ri—§- 0.0803. (2.13) 

The results obtained in this section are summarized in Fig. |2|. The exact leading- 
order results for the gg QQQQ and qq QQQQ processes, obtained with MAD- 
GRAPH/MADEVENT |l|, ||], are shown (histo grams) as a function of y = log{s / ArriQ) . 
The dominance of the i-channel gluon exchange contribution, present only in the case of gg 
initial state, is apparent. The diamonds are obtained by integrating the high-energy limit 
of the matrix element, Eq. ( |2.3D , with the exact phase space; the difference with the exact 
result is fairly small, which implies that, at the dynamical level, the high-energy limit is 
a good approximation. The approximation of the phase space is evidently more drastic. 



and results in the constant (dashed) red line, whose value is taken from Eq. ( 2.13 ). For 
comparison we also show the Born gg, qq — > QQ processes corresponding to the ooo con- 
tribution in Eq. (^). As argued above, and in contrast to the gg — > QQQQ contribution, 
these exhibit a exp(— y) behaviour in the high-energy (large Y) limit. 



-6- 



3. The BFKL Monte Carlo 



As we have seen, in the high-energy hmit ( |2.2[ ) the cross section for the production of four 
heavy quarks is dominated by processes with a gluon exchange in the crossed channel. In 
that hmit, the BFKL formahsm resums the universal leading-logarithmic (LL) corrections, 
of 0(0" log"(s/|t|)), with t defined in Eq. (p.4|). These are obtained in the limit of strong 
rapidity ordering of the emitted gluon radiation, 

yQi - 2/(32 ^ 2/1 ^ 2/2 > • • • » Vn-l > Vn > VQ-s ~ Vq^ , (3.1) 

where we label by 1, . . . , n the emission of n gluons along the BFKL ladder. Because of the 
strong rapidity ordering, the contribution of the gluons to the parton momentum fractions 
( |2.5| ) is subleading, and it is therefore neglected to LL accuracy. The BFKL-resummed 
cross section for the production of four heavy quarks is then given by Eq. (p.6|), where the 
6 function, (5^(qa^ — qfe^)/2, is replaced by the solution of the BFKL equation, 

/(q,^,q,^,Ay) = == ^ e-<^ / dz.e-(-'")^^ , (3.2) 

(27r)2./«2 n? „_ „ ^-00 \Qt. I 



with (p the azimuthal angle between Qa and qb, and uj{i/,n) the eigenvalue of the BFKL 
equation with maximum at ^(0,0) = 41og 2CyiQs/7r. Thus the solution of the BFKL 
equation resums powers of Ay, and rises with Ay as /(qaj^,qfe^. Ay) ~ exp(a;(0, 0)Ay). 

However, in a comparison with experimental data, it must be remembered that the LL 
BFKL resummation makes some approximations which, even though formally subleading, 
can be numerically important: a) the BFKL resummation is performed at fixed coupling 
constant, and thus any variation in the scale at which as is evaluated appears in the next-to- 
leading-logarithmic (NLL) terms; b) because of the strong rapidity ordering any two-parton 
invariant mass is large. Thus there are no collinear divergences in the LL resummation in 
the BFKL ladder; c) finally, energy and longitudinal momentum are not conserved, since 
the momentum fractions x of the incoming partons are reconstructed from the kinematic 
variables of the four heavy quarks only, without including the radiation from the BFKL 
ladder. Therefore, the BFKL theory will severely underestimate the correct value of the 
x's, and thus grossly overestimate the gluon luminosities. In fact, if four heavy quarks -|- n 
gluons are produced, the correct evaluation of the x's yields 

x. = v^^%f^+v^ x. = y "^^-;:''' + vg^, (3.3) 
k ^ k ^ k ^ k ^ 

where p-,_L are the transverse momenta of the gluons produced along the BFKL ladder. 

In the standard (analytic) approach to BFKL, which leads to Eq. (p.2|), it is not possible 
to take the contribution of the BFKL gluon radiation into account in Eq. ( |3.3D . This is 
because in deriving Eq. (^^) one has already integrated over the full rapidity ordered 
phase space for BFKL gluon radiation. To gain information on the BFKL gluon momenta 
we need to unfold the gluon integrations. This approach results in an explicit sum over 
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the number of emitted BFKL gluons, where each term in the sum is an integral over 
the rapidity ordered BFKL gluon phase space. The solution to the BFKL equation can 
then be obtained (numerically) while maintaining information about each emitted gluon 
by evaluating these integrals in a Monte Carlo approach |16, 17|. Besides allowing energy 



and momentum conservation to be observed by including the BFKL gluon contribution 



to Eq. (3.3), this approach also allows subleading effects originating from the running of 
the coupling to be taken into account. The method has recently been generalised to solve 
the BFKL equation at full NLL accuracy |l8| , |l^, although some work remains to be done 
before it can be applied in a phenomenological study like the one presented here. 



The Monte Carlo formulation of Ref. |17] is, in its simplest form, applicable only when 



the transverse momentum of at least one end of the BFKL chain is kept bigger than some 
cut-off |qj| > Pj_ » /u with i E {a, 6}, and fi the resolution scale of the BFKL Monte 
Carlo (see Ref. [|^] for further details). It was demonstrated in Ref. [0] that in the case 
of hadronic dijet production with a minimum Pi = 20 GeV, the residual ^-dependence 
is negligible for fi < 6 GeV. Varying n will shift contributions between different /("^'s 
describing the contribution from different numbers of resolved gluons. 

However, in the current process of 4Q production there is no minimum transverse 
momentum scale at either end of the BFKL chain. To resolve the problem thus faced 
by the BFKL MC formulation we cut out a small region of phase space corresponding to 
p± < 0.05 GeV at one end of the chain. The contribution from this very small region of 
phase space is negligible, but nevertheless this cut-off is sufficient to permit the use of the 
unfolded BFKL formalism. In principle fi could then be chosen arbitrarily small compared 
to the cut-off, but this would result in very slow convergence due to the extremely large 
number of resolved gluons with a transverse momentum above this scale. Instead, fi is 
chosen according to the transverse momentum at one end of the BFKL chain in 5 steps. 
This keeps the average number of resolved gluons under control and thus ensures rapid 
convergence, while maintaining the very weak //-dependence of the overall result. 

In order to demonstrate the behaviour of the BFKL ladder, we consider the production 
of four heavy quarks assuming that all of them are detected. We study the production rate 
as a function of the transverse momentum q^j^ = — Pg^j^ — Pqjx exiting from the impact 
factor I^^{qa). At leading order, the transverse momenta of the two pairs are equal. 



q^^ = qbj^ = qj^. Since we know from Eqs. ( B.13 ) and ( B.14 ) that the scaling of the 
integrated impact factor is g{^) ~ 0{^), with ^ = q^/rnq, power counting from Eq. ( 2.12| ) 
shows that at leading order da/dq^ ~ (^(q^) Q± ~^ 0- When the BFKL gluon radiation 
is included, the production rate is hardened in the infrared and we obtain 

^'^ const. as qa^^ — > , q^^ . (3-4) 



dqa,dqb 



In Fig. ^ we plot the transverse momentum distribution da / dqa^dq^^ evaluated at qa^ = 
qiy^ = . The solid red curve is the four-quark production (|2.6| ) , but with the high-energy 
parton momentum fractions replaced by the exact ones, — > Xa and — > x;,; in this case, 
the two impact factors have equal transverse momenta q^^ = q^j^ . The dashed blue curve 
corresponds to the high-energy limit of leading-order four-quark production (|2.6|) with the 
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1.2 



LO (HE) 




pJGeV] 



Figure 3: The transverse momentum distribution da/dqa^dqi,^ evaluated at dqb^^ = dqt^ = P±- 
The solid red curve corresponds to the high-energy limit of leading-order four 6-quark production, 
with nib = 5 GeV. The dashed blue curve corresponds to adding BFKL evolution to the gluon 
exchanged in the i-channel. 



BFKL ladder included. In this case is no longer restricted to be equal to q;,^ , which 
explains why the BFKL curve is lower than the leading-order one. However, we see that 
the spectrum is relatively harder for p_L — > in the BFKL case. 



4. BFKL signals at the Tevatron and LHC 
4.1 Inclusive heavy-quark production 

In this section we compare the results for the 4(5 channel, obtained with the BFKL MC 
described in the previous section, with those relevant to QQ production, obtained with the 
NLO code of Ref. ||5| and MC@NLO |2^. We consider bottom quark production, with 



mf, = 5 GeV, since 6-quarks are readily identifiable at the Tevatron and LHC. In the case 
of pair production, we need to use a NLO computation in order to explicitly verify that, 
with our chosen set of cuts, large non-BFKL logarithms do not appear in the cross section, 
which is a necessary condition in order to study BFKL signals with the 4Q channel. 
Figure Q shows the integrated cross section 

a{Ay) = / dAy'-—{Ay') (4.1) 
JAy dAy 

as a function of Ay, the rapidity distance between the two tagged quarks (which, for this 
process, are b and b), at Tevatron and at LHC energies. In order to simulate a realistic 
detector coverage, the rapidity of both quarks is required to be less than 2.5, and therefore 
Ay = 5 is the largest accessible rapidity separation. We also consider additional cuts on the 
transverse momenta of the tagged quarks, imposing prp^ 5 > 5 and 10 GeV. The two-loop 
running of the strong coupling as, and the MRST99 package of pdfs has been used, with 
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Figure 4: Integrated cross sections as a function of Ay, at LO (dashed histograms) and NLO 
(solid histograms), with no cut on the transverse momentum p^, and with px > 5 and 10 GeV, at 
the Tevatron (left panel) and LHC (right panel) energies. The code of Ref. j|] has been used. 



factorisation scale set to = (m^^ + m| )/2. ^From the figure we can see that the cuts 
on the transverse momenta largely reduce the impact of radiative corrections. However, 
this information alone is not sufficient to guarantee that non-BFKL logs do not spoil the 
perturbative expansion. In order to investigate this issue, we thus recomputed the cross 



section with MC@NLO [gO|, |2l|], which, by matching the NLO results with the HERWIG ^| 
parton shower, improves the fixed-order result by effectively resumming various classes of 
large logs. In the case in which no pT cuts are applied, the MC@NLO results are basically 
coincident with the NLO ones. However, by imposing prpj^ ^ > 5 GeV, the MC@NLO cross 
section is roughly a factor 1.7 larger than the NLO, in the whole Ay range considered. This 
is due to the fact that the pT cuts render the cross section sensitive to Sudakov effects. 
Although these could be reduced by imposing different pT cuts on the two tagged 6's, it 
is quite problematic to eliminate them completely. Thus, the pure NLO result must be 
regarded, at least for the pT cuts considered here, as a lower bound on the bb inclusive 
cross section. 

In order to be definite, we require Px6b > ^ GeV in what follows. In Fig. [s], we 
plot the integrated cross section for bb production as a function of Ay, at Tevatron and 
at LHC energies. For the sake of comparison, we display again here the middle NLO 
curves of Fig. ^. In addition, the dot-dashed red curve displays the high-energy limit 
contribution of the AQ channel to inclusive 2Q production, where Q = b ox b. The fac- 
torisation and renormalisation scales have been set to = = {jnl^^ + m|^^)/2 and 
/^Ffe ~ ^^\b ~ ("^bsx ~'~ ™'54j^)/^' Thus, the strong coupling a| must be understood here as 
Q;|(/i|„) a|(^|j), with as evolved at two loops, in accordance with the NLO calculation. -t 
The dashed blue curve is the same as the red curve but with the addition of BFKL gluon 



'''We justify the scale choices as follows: in the high-energy limit the impact factors for hb production on 
either side can be viewed as two almost independent scattering centres linked by a gluon exchanged in the 
crossed channel. It therefore makes sense to run the pdfs and as according to the scales set by each impact 
factor. 
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Figure 5: Integrated cross sections as a function of Ay, with px > 5 GeV, at Tevatron (left panel) 
and LHC (right panel) energies. The NLO curves are the same as the middle NLO curves of Fig. 4, 
and are displayed here for the sake of comparison. The dot-dashed red curves are the high-energy 
limit contributions of the AQ channel to inclusive 2Q production; the dashed blue curves are the 
same as the red curves with the addition of BFKL gluon radiation. 



radiation. In the BFKL gluon emission chain, the value of is taken at the b mass. 

Fig. |5| is the central result of this study. It shows that within the rapidity range for 
heavy-quark production accessible to LHC (assumed here to correspond to Ay < 5) the 
4Q channel, even augmented by the BFKL gluon radiation, can never overcome the 2Q 
channel. Thus, it cannot readily be used as a footprint of BFKL radiation. 

The situation could be improved either by imposing the additional requirement that 
the two tagged 6-quarks have the same sign, i.e. bb or bb, or by requiring three or more 
6-quarks to be identified. In the former case, this would reduce the "4(5" curves in Fig. ^ 
by a (combinatoric) factor of two, while almost completely removing the 26 contribution. 
A realistic assessment of how much of the BFKL signal would remain in these cases would 
depend on the efficiencies of multi-6-quark tagging and charge identification (via the sign 
of the lepton in semi-leptonic i?-meson decay, for example), which goes beyond the scope 
of the present study. Another issue that needs to be addressed by a more realistic study is 



the contamination from overlapping events, see for example Ref. |24) 



4.2 Inclusive heavy-quark -|- 1 jet production 

As mentioned in the Introduction and in Section |2|, inclusive QQ+1 jet production is also of 
interest in the high-energy limit, and is in a sense a hybrid of the original Mueller-Navelet 
2 jet and our 4Q processes. In this process, a gluon is exchanged in the t-channel between 
the jet and the QQ pair already at leading-order, which in this case is 0{al). In fact, as 
in Eq. (p.l|), the Ay distribution for QQ + 1 jet production in the high-energy limit can be 
written schematically as 

d(j - ^ °° 

~ 4Y,boM + 4Y,biMsLy + 4Y,b2jas{asLy +•••, (4.2) 

y 7=0 1=0 7=0 
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Figure 6: Inclusive heavy two-quark +1 jet production as a function of the rapidity separation of 
the jet from the average position of the heavy two-quark pair, Ay = \yj — {yi, + yi,)/2\. The soHd 
and red dot-dashed curves correspond to leading-order production, exactly and in the high-energy 
limit respectively. The dashed blue curve is the leading order plus BFKL rcsummation. 



where L = log(s//i^) ~ Ay is a large logarithm, and the quantity //^ is a mass scale 
squared. The first sum in Eq. (4^) is a fixed-order expansion in as starting at 0{ag), and 
collects the contributions which do not feature gluon exchange between the jet and the 
QQ pair. The boo coefficient is the leading-order term (a specimen diagram is depicted in 
Fig. |l|(b)). The second and third sums of Eq. (4.2) collect the contributions which feature 
only gluon exchange in the crossed channel between the jet and the QQ pair, the second 
(third) sum resumming the BFKL (next-to-) leading logarithmic corrections. Fig. |l](c) 
represents the zeroth-order term of the second sum. The bij and b2j coefficients behave 
like l//^tw- We note, however, that in contrast to Eq. (^3), the second and third sums 
of Eq. ( |4.2D start at the same order in Qj as the first sum. Thus one would expect that 
the onset of the dominance of the asymptotic energy region s — > oo occurs more quickly 
in this case than in heavy two-quark production. We analyse this issue by computing the 
coefficients 5oo and bij. 

We consider inclusive heavy two-quark -|- 1 jet production in the high-energy limit. 



The heavy quarks are b quarks for which, following the analysis of Section 4.1, we require 



that \yi,i\ < 2.5 and PTfe,fe > ^ GeV. For the jet, we require the set of cuts \yj\ < 5 and 
PTj > 30 GeV. The factorisation and renormalisation scales are taken as ^^^j = /i 



2 

Ra 



("^bix ~'~ ^^'pb ~ ~ Ptj^ since the impact factors for bb production on one 

side and for jet production on the other can be viewed as two almost independent scattering 
centres linked by a gluon exchanged in the crossed channel. Thus the strong coupling a| 
must be understood here as a|(/x|^) as(/i|j). In Fig. ^, we show the distributions for 
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heavy two-quark + 1 jet as a function of the rapidity separation of the jet from the average 
position of the heavy two-quark pair, Ay = \yj — (yb+Vij/'^l ■ The sohd curve corresponds to 
leading-order production (exact matrix element); the dot-dashed red curve is the leading- 
order production in the high-energy limit approximation and the dashed blue curve is the 
leading order plus BFKL resummation, as given by the Monte Carlo generation of the 
ladder gluons (i.e. with energy- momentum conservation). 

There is evidently a sizeable suppression from the resummation when the BFKL glu- 
ons are radiated off the ladder while conserving energy-momentum, reminiscent of what 
happens in the case of dijet production in the high-energy limit [25|. This is at first sight 
puzzling, because the kinematics of two-quark -|- 1 jet production in the high-energy limit 
resemble more closely the ones of W+ 2 jet production rather than those of dijet production, 
and in W + 2 jet production in the high-energy limit there is no such strong suppression 
when enforcing energy-momentum conservation on the BFKL ladder pi. However, that 
is where the similarity ends: in 1^ -|- 2 jet production, the impact factor for 1^ -|- 1 jet 
production is generated by a quark, while in the present case the impact factor for bb pro- 
duction is generated by a gluon, and therefore the dependence on the pdfs in the two cases 
is completely different.^ 



5. Conclusions 

A definitive test of BFKL physics at hadron colliders is still lacking. A number of processes 
have been suggested, including the standard Mueller-Navelet dijet production, and in this 
paper we have studied a new possibility: four heavy-quark production with a large rapidity 
separation between two of the heavy quarks. The common feature of all these 'BFKL' 
processes is the presence of a t-channel gluon in the scattering amplitude, which gives the 
dominant contribution in the high-energy limit. 

In this work we have focused on the production of b quarks at Tevatron and LHC 
energies. The simplest quantity to measure is the 2b inclusive cross section as a function 
of the rapidity separation Ay. However in this case the 46 process has to compete with 
leading- and next-to-leading-order bb production. Using a set of representative cuts on 
rapidities and transverse momenta, we have shown that in practice the NLO bb contribution 
is dominant over the measurable Ay range, although at the very highest Ay values (~ 5) 
at the LHC energy the bb and 46 contributions are of comparable magnitude. 

We can conclude, therefore, that it will be very difficult to detect any BFKL signal in 
the 26 inclusive distribution. However, a characteristic feature of the 4Q contribution in 
the high-energy limit is that the two heavy quarks separated by a large rapidity distance 
are as likely to have the same as opposite sign. The ability to tag the sign of the 6 quarks 
could therefore be used eliminate the NLO 66 contribution. We note also that in the case 

^In order to rule out other possible explanations, we tried to mimic as much as possible the set-up of 
W + 2 jet in two-quark -I- 1 jet production, namely we eliminated the gluon-gluon sub-process, so as to 
make two-quark + 1 jet production by quark-gluon scattering the dominant process, and we set the &-quark 
mass equal to the W mass. Even with these modifications, we still obtain a BFKL distribution with same 
qualitative features as in Fig. n. 
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of the 45 process, the bottom quantum number is conserved locally in rapidity, i.e. many 
of the events with two detected b quarks with a large rapidity separation could have one 
or two additional b quarks in the detector. To study these possibilities in detail would 
however require detailed knowledge of the detector capability, and is therefore beyond the 
scope of the present work. 

Finally, we also considered the case of QQ + 1 jet production, which is an extension 
of the original dijet case in which one of the far forward/backward jets is replaced by an 
heavy-quark pair. Here there is a i-channel gluon already at leading order and so one 
might expect an earlier onset of the high-energy asymptotic regime. However, because 
the dominant contribution involves gluons in the initial state, there is a severe suppression 
from the pdfs when the additional energy radiated in the BFKL ladder is properly taken 
into account. This means that for this process, the fixed-order pcrturbative contribution 
(i.e. LO or NLO) is likely to be a good approximation to the full cross section over the 
accessible kinematic range. 
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A. Four heavy-quark production 

The partonic cross section for four heavy-quark production is 

da{PaPb PQiPq^PQsPqJ = ^dV^lMp^p^^Q^Q^Q^gf , (A.l) 

with s = (pa + Pb)"^ = XaXjjS and S the squared partonic and hadronic centre-of-mass 
energies respectively, and with four heavy-quark phase space 

= n iS^tfi- (2^)' ^^^P- + f f - PQ^ - PQ^ - PQs - PqJ > (A.2) 
i=i ^"^^^ ^PQi 



with Pq . = y^P^~+^^- The factorisation formula is 

da = ^ dxadxb fa/A{xa, I^V) fb/B{xb, f^p) da , (A.3) 

ab 

where the sum is over the parton species, and fa/A{xa, l^'p) is the pdf of the parton a of 
momentum fraction Xa within hadron A, and similarly for parton b. Parametrizing the 
heavy-quark momenta in terms of the rapidities, 

PQi = (mQi^ cosh yQ, , pg.^ , mg^^ sinh yg. ) , (A.4) 
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we can write the cross section for heavy-quark production as 
da 



E 



x^xi (2i)-i(47r)iS2 

with momentum fractions of the incoming partons given by 



X. 



i=l 



Vs 



Xb 



E 

i=l 



-VQi 



Vs 



(A.5) 



^ (PQix + PQ2X + PQ3X + PQ4x) ' 



(A.6) 



B. Impact factor for gg* QQ 

The impact factor, for QQ production can be obtained by using the squared amphtude 
for gi ^ QQi with i = q,g from Refs. ^. The momenta of the incoming and outgoing 
partons are g{pa) + i{Pb) = Q{pq) + Q{Pq) +iiPb')- In the high-energy hmit, the rapidities 
are strongly ordered while the transverse momenta are of similar size, 



(B.l) 



yQ^yQ^Vi, PQ± - pq± ~ Pij- ■ 

The squared amplitude for gi ^ QQi, summed (averaged) over final (initial) colours and 
helicities, then reduces to 



\M 



gi-*QQi 



4^2 

^ = -^I^^iPa,PQ,PQ;q)PiPb,Pb') , 



q,9: 



(B.2) 



with s = {pa + Pb)'^ the squared centre-of-mass energy, and t = {pb — Pb')^ the momentum 
transfer. The impact factor P{pb,pb') for quark/gluon production, summed (averaged) over 
final (initial) helicities and colours, can be written as ||2^ 

-.2 



9 



(B.3) 



A^2 _ 1 ' - 2Nc 

where Ca = Nc = 3, the index c runs over the colours of the gluon exchanged in the 
crossed channel, and we have used the standard normalization of the SU(A'^c) matrices, 
tr(A'^A'^') = 5^^' /2. The impact factor for gg* QQ, summed (averaged) over final (initial) 
colours and helicities, is then 



,4 xcc' 



.Pa,PQ,PQ;' 



9t 



t' 



i (1 + 2iV2 ^ ^ 



/^2 4aQf. 



t' 



V 2 _ 

QQ 



+ ^ I ^*'aQ ~ ^^'aQ 



+ {Q ^ Q,X ^ x) , 

where we have defined the momentum fraction 



(B.4) 



X 



+ p 



1 — X . 



Q 



(B.5) 
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and the invariants 



m 



+ 



Q± 



t'aQ = (Pa-Pq) 



m; 



2 h_^- 
Q - 



m: 



Q± 



t' 



{Pa 



m: 



2 h_^- 



m 



Q± 



(B.6) 



In the small qj_ limit, the jet opposite to the impact factor for QQ production becomes 



collinear, and the cross section obtained from the squared amplitude (B.2) yields an infrared 
singular real correction. Since the latter may have at most a logarithmic enhancement as 
q_L — > 0, the squared amplitude ( |B.2| ) cannot diverge more rapidly than l/q\. This means 
that in the small q±_ limit, the impact factor must be at least quadratic in q^, I^^ ~ 0{q\). 
Using = — (pq±+Pq^), we see immediately that this is the case. In addition, as — > 
we have an almost on-shell gluon scattering with a gluon, then pg^ — > — Pq_l and averaging 



over the azimuthal angle of qj_, Eq. ( B.4 ) becomes 



lim 



at 



q\ xi {N^ - 1) 



2N^xx 



1 - 2xx 1 



2m^ 



m 



q± XX 



m 



\M 



Q± 
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(B.7) 



where for the invariants in the g g ^ QQ Born amplitude we have used 



^q+Pq± 

X X ' 
X^Q +Pq± 



(B.8) 



2m: 



Q 



t-s. 



Equation (B/7) explicitly shows that in this limit the impact factor is positive definite, and 
that it factorises into the squared amplitude for g g ^ QQ scattering. 



B.l The integrated impact factor for gg* QQ 



Using Eq. (B^) and the invariants ( B.6 ), the integrated impact factor ( 2.11 ) becomes 



2 xcc' 



aid' 



2iV,(iV2 - 1) 
(x2 + x2)k2 



+4xxmQ 



dx 



1 



TT 

+ TZ-^ + 



D1D2 DiD-i D2D3 
1 N}. ^ 



D1D2 DiD3 D2D3 



2Dl 



1 + 



1 N' 



2Dl 



+ 



D2 

-^3 



(B. 



(B.IO) 
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where on the left-hand side we have made exphcit that the impact factor depends only on 
the dimensionless ratio ^ = g^/mg. In Eq. ( B.10| ) the propagators 



Di=mQ+ p\ 



(B.ll) 



have been used. Introducing the Feynman parameter A and performing the integration 
over the transverse momentum gives 



with 



5(0 



1 



2N,{N^ - 1 



1 .1 
dx / dX 
Jo 

-^ + 2ix{2 + 2N^ X [-Ai + X [x"^ + x"^) C] 



(B.12) 



1 + (1-A)AC 



+ 



l + (l-A)x2Ae 



. (B.13) 



Note that as ^ ^ 0, g{0 ~ 0(0, in accordance with Eq. ( |B.7] ). As ^ — > cc, it grows 
logarithmically, g{(,) ~ log(0- The integrals of Eq. ( p.l3| ) can be performed analytically, 
and we obtain 



9(0 



1 



9iv,(iv2 - i)e 

-4N^ {5C - 12) - 3^ + 



(B.14) 



12 [2jV2(g-2)(e + 4)-Oe-l) 
v/Oe + 4) 



■ tanh 



C + 4 
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